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Abstract 

We consider a different power counting in potential NRQCD by incorporating 
the static potential exactly in the leading order Hamiltonian. We compute the 
leading relativistic corrections to the inclusive electromagnetic decay ratios in 
this new scheme. The effect of this new power counting is found to be large (even 
for top). We produce an updated value for the decay to two photons. This 
scheme also brings consistency between the weak coupling computation and the 
experimental value of the charmonium decay ratio. 



1 Introduction 



The determination of heavy quarkonium properties from QCD has always been 
a major objective in high energy physics. In this respect, the development of 
effective field theories (EFT) directly derived from QCD like NRQCD or 
pNRQCD |2j (for a review see Ref. [3]) has opened the door to model inde- 
pendent determinations of heavy quarkonium properties. Instrumental in this 
development is the fact that heavy quarkonium systems can be considered to 
be non-relativistic (NR). They are then characterized by, at least, three widely 
separated scales: hard (the mass m, of the heavy quarks), soft (the relative mo- 
mentum IpI ~ mv <^ m, of the heavy-quark-antiquark pair in the center of mass 
frame), and ultrasoft (the typical kinetic energy E ~ mv"^ of the heavy quark in 
the bound state system). 

In this paper we focus on pNRQCD. This EFT takes full advantage of the 
hierarchy of scales that appear in the system, 

m ^ mv ^ mf ^ ■ ■ ■ , (1) 

and makes a systematic and natural connection between quantum field theory 
and the Schrodinger equation. Schematically the EFT takes the form 



(^i^o-^- v;(°)(ol <f(r) = 

\ m J 



> pNRQCD 



-|- corrections to the potential 

+ interactions with other low-energy degrees of freedom 

where Vj'^\r) is the static potential and $(r) is the Q-Q wave function. 

A major issue to be settled is to decide upon the precise form of Vj'^^r)^ in 
particular whether one works in the weak or strong coupling regime and how to 
treat subleading terms. In the strict weak coupling regime one could approxi- 
mate the static potential by the Coulomb potential Vj^\r) ~ Vc = —CpOL^jr 
and include higher-order terms perturbatively. There seems to be growing con- 
sensus that the weak coupling regime appears to work properly for t-t production 
near threshold, the bottomonium ground state mass, and bottomonium sum rules 
(for a recent discussion on this issue see ^). One would then expect that other 
properties of the bottomonium ground state like the hyperfine splitting or elec- 
tromagnetic decay widths could be described as well by the weak coupling version 
of pNRQCD. However, in this case the situation is not that clear. There has been 
a precise determination of the bottomonium ground state hyperfine splitting us- 
ing the renormalization group in pNRQCD [5]. Nevertheless, the predicted value 
does not agree well with the recently obtained experimental number |6l [7] . There- 
fore the situation remains unsettled. For the inclusive electromagnetic decays the 
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convergence is not very good [8]. Even for top, higher-order corrections to the 
normahzation appear to be sizable |8l [9l [101 E] • 

In principle, the novel feature of these observables (maybe more so for the de- 
cays) compared to the heavy quarkonium ground state mass is a bigger sensitivity 
to the value of the wave function at the origin and to its relativistic corrections. 
Note that in this case the relativistic corrections are divergent and their diver- 
gences have to be absorbed by the matching coefficients of the effective theory: 
potentials and current matching coefficients. If one considers the decay ratio, the 
dependence on the wave function associated to the static potential drops out and 
only the relativistic correction survives. This makes the decay ratio the cleanest 
possible place on which to quantify the importance of the relativistic corrections 
to the wave function. 

In Ref. [12] the decay ratios have been computed with NNLL accuracy, ac- 
counting for the resummation of logarithms. The scale dependence has greatly 
improved over fixed-order computations and the result is much more stable. The 
convergence could be classified as good for the top case, reasonable for the bot- 
tom, and not good for the charm, although in all three cases the scale dependence 
of the theoretical result was quite small. For the case of the charm there is ex- 
perimental data available, and the agreement with experiment deteriorates when 
higher order corrections are introduced. On the other hand there exists a nice 
analysis for charmonium in Ref. |T3], where they consider a potential model (a 
Cornell-like one, yet compatible with perturbation theory at short distances, since 
it is coulomb-like in this regime) for the bound state dynamics, but a tree-level 
perturbative potential for the spin-dependence. They also correctly performed 
the matching in the ultraviolet with QCD along the lines of what would be pN- 
RQCD in the strong coupling regimo Their net result was that they were able 
to obtain consistency with experiment albeit with large errors. Unfortunately, 
this result suffers from model dependence. In particular, since a perturbative po- 
tential has been used for the spin-dependent potential, it would have been more 
consistent to treat the static potential also in a perturbative approach. In this 
respect, it has been shown in Refs. [HI [151 HSl HZ] that, once the renormalon 
cancellation is taken into account, the inclusion of perturbative corrections to the 
static potential leads to a convergent series and that this series gets closer to the 
lattice values in the quenched approximation up to scales of around 1 GeV. It is 
then natural to ask whether the inclusion of these effects may lead to a better 
agreement in the case of charmonium and for sizable corrections in the case of 
bottomonium and t-t production near threshold. Note that in this comparison 
between lattice and perturbation theory one has to go to high orders to get good 
agreement. Therefore, a computation of the relativistic correction based on the 

^Actually the whole computation would fit into the strong coupling regime of pNRQCD 
except for the fact that the spin-dependent potential is computed in perturbation theory. 
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leading order expression for the static potential, i.e. the Coulomb potential, as 
the one used in an strict NNLL computation, may lead to large corrections, since 
these corrections, as well as the wave function at the origin, could be particularly 
sensitive to the shape of the potential. 

Therefore, in this paper we reorganize the perturbative expansion and con- 
sider the static potential exactly, whereas we treat the the relativistic terms as 
corrections. By doing so we expect to have an effect similar to the one observed 
in Ref. [13]. Including also the renormalization group improved expressions, we 
expect to obtain results with only a modest scale dependence. The explicit com- 
putation will confirm to a large extent these expectations. We will be able to 
give an updated prediction for the decay of the rjb to two photons and obtain a 
result for the charm decay ratio compatible with experiment (though in this last 
case with rather large errors). Note that our computation is completely based 
on a weak coupling analysis derived from QCD and no non-perturbative input is 
introduced. 



2 Decay ratio 

The one-photon mediated processes are induced by the electromagnetic current 
jfj,, which has the following decomposition in terms of operators constructed from 
the non-relativistic quark and anti-quark two-component Pauli spinors and x 

m- 

j = cM^^ax + + • • ■ , (2) 

where is the renormalization scale, D is the covariant derivative, cr is the 
Pauli matrix, and the ellipsis stands for operators of higher mass dimension. The 
Wilson coefficients c„ and represent the contributions from the hard modes 
and may be evaluated as a series in in full QCD for free on-shell on-threshold 
external (anti) quark fields. We define it through 

= Ef^y^^H/^) , c(°)=i, (3) 



i=0 



IT 



and similarly for other coefficients. The coefficients c[,^^ and c[,^^ have been com- 
puted in Refs. [19] and [20l El] respectively. 

The operator responsible for the two-photon S'-wave processes in the non- 
relativistic limit is generated by the expansion of the product of two electromag- 
netic currents and has the following representation [18] 



O,, = c,,{f,Wx + ^^^^D\ + . . . , (4) 



6m^ 
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which reduces to the pseudo-scalar current in the non-relativistic hmit. The 
coefficients ci^^ and c^^ have been computed in Refs. [22] and [23] (in semi- 
numerical form) respectively. 

Let us define the spin ratio for the production and annihilation of heavy 
quarkonium Q as 

' a(77 ^ Qin'So)) T{Q{n^So) ^ 77) ' ^ ^ 

The effective theory expression for the spin ratio reads 

where is the quark electric charge, Cs{fi) = Ct,(/i)/c^^(yu), ^/'^"'^^(r) are the spin 
triplet (vector) and spin singlet (pseudo-scalar) quarkonium wave functions with 
principal quantum number n. The wave functions describe the dynamics of the 
non-relativistic bound state and can be computed within pNRQCD. The latter 
is the Schrodinger-like effective theory of potential (anti) quarks whose energies 
scale like m^f ^ and three-momenta scale like rUgV, and their multipole interaction 
to the ultrasoft gluons [211 [251 [261 [27]. The contributions of hard and soft modes 
in pNRQCD are represented by the perturbative and relativistic corrections to 
the effective Hamiltonian, which is systematically evaluated order by order in a.^ 
and V around the leading order (LO) Coulomb approximation. 



3 pNRQCD framework 

As we have mentioned before, the framework we use to compute the decay ratio, 
and more specifically the wave function, is pNRQCD. For the purposes of our 
paper the full setup of pNRQCD is not needed. We will only need the static 
potential, Vj^\r), and the spin- dependent potential Vg2\{r). Furthermore, we 
will reorganize the perturbative expansion. The static potential will be treated 
exactly by including it in the leading-order Hamiltonian 

/7W^-^ + yi°)(r), (7) 

where rrir = mim2/{mi + m2)- On the other hand, the spin-dependent potential 
(in D = 1 + (i = 4 — 2e dimensions) 

Aif = = --r^^ [SI, Si][S^, S^#W(r) (8) 

mim2 0,17111712 
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is considered to be a perturbation to the result obtained with H^'^^ Therefore, 
we distinguish between an expansion in v and Og- v has an expansion in ag itself 
but this expansion does not converge quickly for these relativistic corrections. 
This remains so even after the inclusion of the renormalon cancellation, which 
has only a minor impact on the determination of the wave function. This is the 
reason we choose to take the static potential exactly. 

As mentioned in the introduction there are different options on how precisely 
to treat Vj'^^ and we will discuss in Section |6] the various options we consider. 
Roughly speaking we will take the static potential up to NNNLO including also 
the leading ultrasoft corrections. We will also need to define a scheme of renor- 
malon subtraction. Therefore, the general form of the static potential will be 

V}'\r) = VsD{r) + 2 6mx, (9) 

where 6mx represents a residual mass that encodes the pole mass renormalon 
contribution and X stands for the specific renormalon subtraction scheme. We 
will show some specific examples in Section 16.31 In Eq. (Q, Vsd is the short 
distance behavior of the static potential, which is independent of the scheme 
for renormalon subtraction (even if we use a non-perturbative potential). In 
momentum space it reads 

~ ~ A-kCf a, .(0) (g) 
lim = Vso{q) = (10) 

with 5,, (())(?) ~ Ois^jj) (for the precise relation see Eq. ( 15T|) ). where a;s(/i) is the 
QCD coupling constant in the MS-scheme. 

For the spin-dependent potential in momentum space we have 
V^l\fi) = [Sl,Si][S-,S^,] 

1(2) 



+H2-3^ ^''^ 

where 81^2 is the spin operator for heavy quark and anti-quark, respectively and 
D^g2 gilJ') = Cisin) + . . .. In the second line in Eq.l ITT]) the spin projection has been 
done, resulting in S*^ = and 2 for spin-singlet and spin-triplet states, respectively 
(this expression actually corresponds to the regularization prescription of [13] 
for the spin- zero states). We have to keep the term of C(e) because the spin- 
dependent potential generates 1/e divergences. The renormalization procedure 
for these 1/e will be discussed in the next section. 
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4 Wave function ratio 



We now turn to the computation of 



= Pn.(/i) = l + 5p„(/i), (12) 



K(0)P 

Applying Rayleigh-Schrodinger perturbation theory to the problem we obtain 

mim2 

where '^^'^^(0) is the wave function for the LO Hamiltonian H^'^'^ and G{E^^) is 
the reduced Green function aX E = which is defined by 



rJm = v^f (0) - G(4°))^^ v^f (0) + o my , (13) 



The prime indicates that the sum does not include the state n and 

G{E) = GiO, 0; E) ^ lim G{r, r; E) = lim(r| ^^^^ _ |r) (15) 

is the zero-distance limit of the Green function G{r,r'; E), which is the solution 
of the Schrodinger equation 

V2 -I 

+ Vj^\r)-E Gir,r';E) = 5{t-t'). (16) 

2mr J 

The short distance behavior of the static potential Vj'^^r) ~ 1/r makes G{E) 
and, therefore, 6pn divergent. Thus we will need to regularize the Green function 
and we will deal with two different ways to do this: dimensional regularization 
and finite-r regularization. We start by considering the former and will come 
back to finite-r regularization in the next section. 

The divergences in Spn are cancelled by divergences in the Wilson coefficient 
Cg(p). Since the latter have been computed in dimensional regularization we will 
need G{E) in dimensional regularization as well. We denote the corresponding 
bare and reduced Green functions by G^^\E) = G^^\0,0;E) and G^^\E^^')) 
respectively. We remark that the LO wave functions (corresponding to H^^^) are 
finite, thus IV-I^^^^HO)!' = IV'^^^^'HO)!' = \¥nKO)\^- 

Using Eqs. f lTT]) - f|T3l) . the bare expression of Spn{p) in dimensional regular- 
ization can be written as 



6 



In order to obtain the MS-renormahzed expression of 5p„, we need to identify 
the divergences of G^^\E^^^). They are the same as those of G^^\E), are inde- 
pendent of E, and can be computed order by order in perturbation theory, since 
they are related to the short distance behavior of the Green function. We thus 
parameterize the divergent and finite terms of G^^\E) and G^^\eI^^) as 

G(^)(i5(o)) 



rrir 

rrir 
2^ 



(0)1 



(19) 



where B^^))(E; jj) and I3^^o^(E^^; jj) are finite in 4 dimensions, but contain terms 

to all orders in Ug/v^ since the bound-state dynamics needs all order resummation 
in ttg. As will be shown, Eq. (15^ . the ultraviolet divergent part can be expressed 
in terms of the (dimensionfull) bare coupling g"^ = Arras /^^"^ as 



Svre 



An 



-2e 



(20) 



will be removed by renormalization. This has to be done consistently with 
the calculation of other parts order by order in the expansion in (in our case 
MS). The divergences are then absorbed in Cg and we can write 



B% (E(°); /i) + ^mrCpa, + 0{al 



.(21) 



3mim2 

This will have to be combined with the MS subtracted matching coefficient Cg(/i) 
in Eq. ^ to obtain the decay ratio. 



5 Green Function in position space 

The main goal of the present paper is to compute &^\E'^^) or, equivalently, 
I3^^o){E^^; /i), with the effect of the static potential included exactly. This calls 
for a numerical evaluation of the Green function rather than pursuing an analytic 
approach. Numerical calculations are most conveniently performed in coordinate 
space. It is here where finite-r regularization comes into play. In Section 15.11 we 
will discuss this regularization and in Section 15.21 we show how to convert the 
Green function obtained in finite-r regularization by matching into the one in 
dimensional regularization. 



5.1 Regularization of the Green function in position space 

The zero-distance Green function with finite-r regularization is simply defined 
as G^'^\E) = G{ro,ro; E), where tq <^ l/{mas)- In order to compute it, we 
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first have to describe how to obtain a numerical solution for the Green function 
G{r,r'; E) in general, given the static potential V^'^\r). Actually the whole pro- 
cedure holds valid for a generic potential (not unbounded from below at long 
distances) that has the correct, perturbative, short distance limitH. According to 
Eq. dH]) renormalon associated affects are power suppressed. Therefore, they will 
not affect properties associated to the r — )■ limit of the potential. 

For the class of potentials described above, the Green function G{r, r'; E) can 
be constructed from the two independent solutions u<(r),u>(r) of the homoge- 
neous Schrodinger equation (our approach follows Ref. [28], see also Ref. [29] ) 



u{r) = 0. (22) 



Here u{r) represents 'U<(r) or ■u>(r'), which are the solutions to Eq. (122|) that 
are regular for r — > and r' — )■ oo respectively. The angular-momentum term is 
dropped in the Schrodinger equation assuming S-wave contribution because the 
limit r, r' — )■ is taken later. The Green function is written as 
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G{ry-E)= r^) for r < r'. (23) 

\2tx J r r' 

The numerical solution at finite r is obtained by solving the Schrodinger 
equation with boundary conditions at short distances. To this end we prepare 
two independent solutions UQ{r) and Ui{r) that are defined by the following initial 
conditions: For 'Ui(r), which we will call the regular solution we set 

Mi(0) = and u'^{Q) = 1 (24) 

so that 

ui(r) =r + 0(r^). (25) 

This completely fixes Ui{r). 

For the non-regular solution, Uo{r), we can not work this way. Whereas we can 
still take Uo{0) = 1, we can not define Uq{0), as it becomes singular. Therefore, 
we first define u'q for small values of r in the following way 

u'oir) = Co{r,) + 2mr r dr'VsD{r') + 0{r), (26) 

Jtc 

where Co{rc) is an integration constant. Note that Tc > acts as a cutoff to 
avoid the denominator-zero of Vsoif')- The total solution then reads (at short 
distances) 



uo(r 



1 + Co(r,) r + 2m, dr' f dr'Vsnir") + C»(r^) . (27) 

Jo Jtc 



^This opens the possibility of using the same formalisni for pNRQCD in the strong couphng 
regime but then we should also consider a non-perturbative potential in Eq. ([5]), albeit with 
the correct short distance limit. 
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This expression can be rewritten as 

yVsD{r') 



1 + Co(r,) r + 2m,r| j'^ dr'Vsoir') - dr'—^^^^ + ^(r^). (28) 



The derivatives u'q ^(r) and UQ^i{r) at small r are used as boundary conditions to 
solve differential equations by, for instance, the Runge-Kutta method. For later 
convenience we take = l/{^e""^) and fix 

Co(r,) = r dr' r' dr"VsD{r") = 2m, P dr' . (29) 

Tc Jo Jrc Jo Tc 

With this choice the 0{r) term of Eq. ( |28l) for uq is a function of In (yue'^'^r) with 
no log-independent terms 

oo 

Uo{r) = 1 + 2m,r ^ Vn ln"(/ie^^r) + ^(r^) . (30) 

n=l 

The coefficients f„, which can be written as an expansion in powers of as(/u), only 
depend on the coefficients a„ of Vsd (see Eq. flS2]) ). i.e. only on the pure short 
distance behavior of the static potential. This choice will turn out to be very 
convenient for the conversion to dimensional regularization, but the final result 
for G{r, r; E) does not depend on this choice. 

From the two solutions Uoir) and Mi(r) we can construct M>(r) and M<(r) as 
follows: First the solution at short distance u^{r) is identified as 

M<(r) = Mi(r), (31) 

because limr_j.o 'Ui(^) = 0. The other solution which satisfies \mir^ooU>{r) = is 
given by 

u>{r) = Mo(r) + 5%(E)Mi(r), (32) 
B'-^lm = -Jim {«o(r)M(r)}. (33) 

From the boundary conditions of mq and Ui it follows that we can mix a Ui- 
component into uo{r). However, the precise choice of uo{r) does not affect M>(r) 
because of the invariance under uo{r) — )■ uo{r) + K,ui{r) with k being an arbitrary 
constant. The zero-distance Green function with finite-r regularization is then 
obtained as 



G^^\E) = ^[AM(ro;/.) + 5^i,(E;/.) 



(34) 



AM(ro;/i) = = --2m,CFa,ln(/ie^^ro) + C(a,'), (35) 
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where the last equahty is a good approximation for /xe'^^ro ~ 1. A^^\ro; jj) 
encodes the divergence of G^'^\E) and plays the role of the 1/e pole of G^^\E). 
It is energy independent because it is related to the overall divergence of the Green 
function. Nevertheless, note that according to Eq. (IMl) we define B^^jff^(E; n) by 

subtracting exactly UQ{ro)/ro, which, depending on the potential, will include 
terms with arbitrary powers of a^- B^^\E; fj,) is computed numerically by solving 
the Schrodinger equation and is independent of the regulator tq. Note however 
that it is scheme dependent, i.e. it depends on the specific condition we use 
for Uo(ro). This dependence cancels between A^^^ and B^''Jo^ such that G^'^^{E) 

is independent of the specific choice for Mo(ro). In analogy to Eq. (fT9|) we also 
define ^ 

(36) 



2n 



(r) 



V 



(r) 

Finally we remark that -B|^(o) is independent of the renormalon subtraction 

scheme used, since G^^\eI^^) and A^'''\ro] jj) are; the latter by the definition used 
in this paper. 



5.2 Conversion to the MS scheme 

Once we have the zero-distance Green function G^'^\E) = G{rQ,rQ; E), where 
tq ^ l/{mas), or more precisely -8^*^(0), we have to convert the result by match- 

ing into the one in dimensional regularization -B|^(o), in order to be able to use 
Eq. fl2T]) . We define the difference 

= 5f ) {E; /i) - B^^j,, {E; /i) = 5^, ^) - b'^I /.) . (37) 

This difference between the schemes can be accounted for by a finite [vq and e 
independent) constant. We also use the fact that the ultraviolet divergent term 
of the Green function is energy independent. This means that is energy 
independent and its perturbative expansion is short distance dominated and can 
be computed order by order in as- 

Note that c^^ does not depend on the long distance behavior of Vj^\ only 
on its short distance behavior, which is universal and dictated by perturbation 
theory, i.e. by Eq. ffTOj) . In particular, the result is independent of the pole 
mass renormalon. Therefore, the value obtained for holds true for a general 
potential (not unbounded from below at long distances) that has the correct short 
distance limit. 

Considering the lowest order approximation of the static potential, the Coulomb 
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potential 



as(/i) 



(38) 



the exact solution for this potential, the Coulomb Green function, is known in di- 



mensional regularization and can be expressed in terms of A = Cp a^j \J —2E /rrir 
as 



47r2 



47i-e-7i3 



According to Eq. ( 1181) we thus get 



rMS 



1 1, (-Sm.,E\ 1 



(39) 



(40) 



Turning to finite-r regularization the Coulomb Green function reads 



TT 



2 rrirCp asTo 



— In (fie^^ro) 



2A 2 



-8mrE" 
1^. 



+ 1 - 7s - V^(l - A) 



whereas Uo('"o) for the Coulomb case is given by 

Mo(ro) = 1 - 2m,. tq In (/ie^-^ro) . 



(41) 



(42) 



The expression stops at 0{as). In the Coulomb approximation there are no Ola"^) 
terms in Uo('"o) and, therefore, in Ay^. Using Eq. fl5^ we then find 



B^^^{E) = 2mrCpas 



1 l^^f-SmrE 



2A 



+ i--fE-ij{i- X) 



(43) 



Note that in an strict NNLO or NNLL computation of the decay ratio this would 
be the only term that should be considered. 

Thus we compute c^^ in an expansion in and obtaiijfl 



MS 



(45) 



This constant can also be obtained from the difference between dimensional- 
and r- regularized computations at finite order in as- At the lowest order it 



^If the constant e^^ were not introduced in Eq. (|35p . would read 



1 



(44) 
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corresponds to the computation of one- and two-loop contributions to the Green 
function in both schemes, by considering the difference of their renormahzed 
pieces. We have checked in an exphcit calculation that Eq. ( I45ll is reproduced by 
the difference of two-loop contributions. 

In order to obtain the C(as) corrections to cj^^ one has to include the 0{a'^) 
corrections to the static potential and compute the associated corrections to the 
Green function in both schemes. In principle, this is possible and partial results 
can be found in the literature. Nevertheless, this would go beyond the aim of 
this work, since it would produce corrections that are anyway unmatched by the 
precision of the hard matching coefficient. 

Finally, for a general potential with the right short distance structure, we can 
combine Eq. f l2T|) with Eqs. f l37|) and fH5|) and write 

^P^^^^ = + Im^Ca. + Oia!)) . (46) 

Once we know the MS expression we can also write 5p„ in different schemes. For 
instance, in the "hard- matching" scheme used in Ref. [12] we have 

^^""'^^^ = -I^^S^^^(^) - 2m.C^a. + Oia',)) , (47) 

which will be relevant afterwards. 

^ (r) 

These results enable us to compute the decay ratio in terms of -B^/o), whose 
determination will be discussed in the next section. 



6 Determination of B L 

In this section we determine B^^Jo^{E^^) in several approximation schemes for 

V^^\ We have already mentioned that our idea is to treat the static potential 
exactly, yet we only know its expression up to three loops. There is some freedom 
on how this truncation is performed. This produces a class of potentials to study, 
which introduces some scheme and scale uncertainties. As we have stressed in 
previous sections, the analysis applies to any arbitrary potential with the cor- 
rect short distance behavior and not unbounded from below at long distances. 
Therefore, in what follows we will consider different approximations to the static 
potential. One quality that they have in common is the renormalon cancella- 
tion. We have to preserve renormalon cancellation between the static potential 
and the pole mass of the heavy quark. At the same time we will be forced to 
consider the resummation of logarithms to reproduce the correct behavior of the 
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potential at short distances. Thus we will have to devise schemes where both 
the resummation of the logarithms and the renormalon cancellation is achieved 
order by order in the perturbative expansion. We illustrate this discussion in 
the following sections, where we show the determination of B^^j^-^ using either 

the Coulomb potential, the static potential at different orders in as{fi), and the 
static potential at different orders in ag(l/r). In this last case we will use dif- 
ferent schemes with renormalon cancellation. The dependence on the scheme of 
renormalon subtraction (potential) may give an estimate of the error, since it is 
also a measure of the dependence on the long distance behavior of the potential. 

Finally, let us note as well that, in order for our computation to make sense, 
the successive approximations to the static potential should be convergent (or at 
least small) themselves. We will check this convergence in this section. 



6.1 Coulomb potential 

If we approximate the static potential by the Coulomb potential Vc we can get an 
analytic solution for by directly working in dimensional regularization. Ex- 
panding G[^\E) as given in Eq. (!39|) around its poles at E^^'' = — mrCi?ag/(2n^) 
we can write 



Gi-\E) = + Gf\Ei^^) + 0{E - E^y) (48) 



n(Ef - E) 

with 



1 1 



■ (49) 



Comparing to Eq. (fT9l) we obtain 



B^^iEl^y) = 2m,C,a, ( -\ In JHh^ + ^ _ + 1 _ ^^n) ] (50) 



and thus (/i) in the Coulomb approximation directly from Eq. fl2T]) . 



6.2 Fixed order V^^^: as{fis) expansion 

The standard way to go beyond the Coulomb potential approximation for the 
static potential is to make an expansion in as(/is). Thus we write 

~ 47rC^a,(/i,) f ^ f as{fis) Y^ , A .... 

ysDKq) = — [~^^ ) ^' ) ■ ^ ' 
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This expanded version of the static potential is often used in quarkonium phe- 
nomenology to respect rigorous expansion according to non-relativistic power 
countingj. In position space we have 

ii^V}«Hr) = Vso(r) = -^£f^{l + £(^^fi^)".„(^.;.)} (52) 



In practice we will take the static potential up to NNNLO, i.e. up to 
including also the leading ultrasoft corrections. This means we take into account 
the first three terms of this expansion with coefficients 

ai{fis; r) = ai + 2/3o In (yUsC^-^r) , 

a2(/i.; r) = as + -jt^o + ( 4ai/3o + 2/3i) In (/i.e^^r) + ifii^ In^ {fisC^^r) , 

571-2 

asifJ's; r) = 03 + ai(3^7r'^ + — /^oA + ISCs/^o^ 

+ (l2ai/3o' + 10/3o/3i) In' (/x.e^^r) + 8(3^ In' (/x.e^^r) 

+ 5a^"(/i„/i„,), (53) 

Explicit expression for ai{fis] r) can be found in the literature [301 EB [321 [331 [SI 
[35l[36]- For the ultrasoft corrections to the static potential we take 

^ar(/^.,M^yC'l7r'ln(^) • (54) 

We will not consider the renormalization group improved ultrasoft contribution 
in this paper as its numerical impact is small. The potential is shown in Figure [1] 
(dashed lines) for /i^ = 2 GeV and the number of light flavors set to Ni = 4. It 
is clear that for small r, depicted in the inset of Figure [H there are serious issues 
regarding the convergence. The potential changes drastically in going from LO 
to NLO to NNLO etc. This behavior occurs for the typical values of fig and A*"; 
that apply for the charm and bottom case. As one increases the value of fis, one 
has to go to shorter distance to see this effect, as it would happen for top. 

Ignoring this problem for the moment and working with the Fixed Order (FO) 
static potential we can obtain uo{ro) and, therefore, A^^lo^{ro) as an expansion in 
as as well. We find 

^So) ('^o) = — - 2m,. ttsifJ.) Cf v{lo) , 



''in the most rigorous fixed order coraputation only the Coulomb part of the static potential is 
treated exactly and as corrections corresponding to the second and remaining terms in Ea.((CT|) 
are treated iteratively order by order by insertion. 



14 



0.0 




1 2 3 4 5 

r [GeV"'] 



Figure 1: The FO (dashed) and RGI (sohd) static potential VsD{r) according to 
Eq. (152]) and Eq. respectively. We take /i^ = 2 GeV, A^; = 4 and //^ = 2 GeV. 
The potential is shown as a function of r at LO (yellow), NLO (green), NNLO 
(blue) and NNNLO (red) with the small r region shown in the inset. The shaded 
area in blue indicates the short distance regime < r < l/m^,. 



'{lo) = E^n(/o) ( 
i=0 \ 



An ) ' 

where lo = In {/i e'^^ro) and the expansion coefficients are given by 
^^o(^o) = ^0, 

Vi{lo) = {ai-2f3o)lo + Pol^ , 

^^2(/o) = - 4ai/3o + 8/3o' + y/3o' - 2/3i) /o 

vsik) = (^as + - 6a2/3o + 24ai/3o' + ai/3o'7r2 - (48 + 27r2) /^ff - 4ai/3i 



(55) 



+ ( 20 + — 1 /3o/3i - 2/32 + 16/3o^C3 - ^''^a ) 



2 

2 




+ (^3a2/3o - 12a^(3^ + (24 + vr^) (3;^ + 2a,(3i - lOM, +(32 + -^cf^ l, 
+ (Aa,P,' - 8/33 + / 3 ^ 2/34 / 4. (56) 
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The /i dependence appearing in Eq. (155|) enters through Eq. (129!) and should be 
cancelled in Eq. ([6]). Even though the exact expression for the static potential is 
scale independent, working at a finite order in as (/is) there is some residual /i^ 
dependence. 

The computation of -B^(o) is done numerically along the lines of Section |5l 

We use the input values ml,ps(2 GeV) = 4.515 GeV [37] and mc,ps(0.7 GeV) = 
1.50 GeV ^38j for bottom and charm quarks, respectively. They can be trans- 
lated into scale-invariant MS-mass of nib = 4.19 GeV and rric = 1.25 GeV. The 
strong coupling ai"'-'^~^\Mz) = 0.118 is used as an input evolved down to low 
energy scale using 4-loop running formulae. For the top quark mass we use 
^t,Fs{'20 GeV) = 173 GeV for illustration. The scale fius needed for the leading 
ultrasoft contribution is set to fius = 0.7 GeV for charm, /i^^ = 1 GeV for bottom 
and fius = 10 GeV for top. 

In Figure |2] we show the results for charm, bottom and top (dashed lines) as a 
function of the scale fig for fixed fi. For illustration we have chosen = 1.5 GeV 
for charm, 2 GeV for bottom and 20 GeV for top. Note that, ideally, the result 
should be independent of fig, as it reflects a dependence on the long distance 
behavior of the potential. For charm and bottom we see problems of convergence, 
in particular for small values of fig- This is due to the behavior of the potential 
at short distances, which we have already illustrated in Figure [H For top the 
situation is much better. Note that the LO curve corresponds to the Coulomb 
potential. In all three cases we observe a significant gap between the Coulomb 
solution and the higher order corrections (for the range of fig for which the result 
can be trusted). 

Before we address the problem of the bad convergence, let us remark that 
expanding the potential in as{fi), the pole mass renormalon enters as an r- 
independent constant in the potential. This constant cancels in the evaluation of 
which is independent of the overall normalization of the potential. Thus, 
in this evaluation the dependence will only enter in the values of the mass used. 
The error associated to this uncertainty is beyond our accuracy. 

6.3 RG-Improved Vj^^: a^il/r) expansion 

In the previous subsection we have seen that the convergence for -B|J{o) is very 
unsatisfactory if we use Eq. f l52p . Surprisingly the problem comes from short 
and not long distances. The solution is to absorb the large logarithms into the 
running coupling. However, this has to be done carefully in order not to destroy 
the renormalon cancellation achieved order by order in as- More specifically, we 
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Figure 2: -B|^(o) as a function of /is at LO (yellow), NLO (green), NNLO(blue) and 

NNNLO(red) with = 1.5 GeV for charm, 2 GeV for bottom, and 20 GeV for 
top. Dashed lines are obtained using Eq. (l52l) . solid lines are obtained using 
Eq. ([59]). 
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consider different approximations to the static potential behaving for r — > as 

- + E a„(l/r; r) (^)" } . (57) 

and yet achieving renormalon cancellation order by order in a"(l/r). This will 
give us an estimate of the dependence of the result on the long distance behavior 
of the potential. We will generically name this class of potentials renormalization 
group improved (RGI) and denote them by LO, NLO, ... according to the power 
of as{l/r) at which we stop the perturbative expansion in Eq. (I37|) . 

One possibility that fulfills all these requirements is the PS scheme [39] with 
the following modification: 

Vpsir) = Vsnir, fir) + 2 6mps (58) 

with 

VsD{r,fir) = / 7^e^^-\/<,^|,=,,(g) + / J^e^^'^Vsn\,=M (59) 

Jq<fj.r [2ny Jq>tir {ZTC)'^ 

Thus we introduce a factorization scale fir- For q < fir we expand Vsoil) in 
as{fis), as in the previous subsection. For q > fir however, we use the running 
coupling in Vsoil)- As can be seen in Figured], the RGI potential (sohd lines) 
shows good convergence for all values of r. We have checked that the results 
for -B'''^(o) are not sensitive to the precise value of the factorization scale, as long 

as fXr is large enough. This definition has the advantage that the renormalon 
contribution is r independent and achieves the resummation of logarithms. The 
fact that the renormalon cancellation is r independent makes it possible to work 
also with 6mps = in Eq. fl58|) . as far as the determination of B^^Jo){E^^) is 
concerned. 

The results for B^^}o)j using Eq. fl59l) rather than Eq. fl52|) are depicted as 
solid lines in Figure [2l We have taken /x^ = 1 GeV, /i = 1.5 GeV for charm, 
fir = fi = 2 GeV for bottom, and fir = fJ' = 20 GeV for top. As can be seen, the 
resummation of logarithms results in a dramatic improvement in the charm and 
bottom case, and is also quite significant in the top case. In all three cases, the 
RG result is nearly independent of fig. This signals a weak dependence on the long 
distance tail of the potential. This is to be contrasted with the results obtained 
using Eq. fl52l) . which are completely unreliable unless unnaturally large values 
for fig are used (especially for charm). The RGI curves show a good convergent 
pattern for top, and also a reasonable convergence in the case of bottom. Even 
for charm we see signs of convergence, albeit marginal. In particular, in this case, 
and to a lesser extent in the case of bottom, the splitting between the NNLO 
and NNNLO curves is not much smaller than the splitting between the NLO 
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and NNLO curves. Note, though, that at NNNLO the potential starts to be 
sensitive to ultrasoft physics, which we do not include in our analysis. In this 
respect the NNNLO curves are to be considered incomplete (though the explicit 
dependence on the ultrasoft factorization scale is small). Moreover, at some point 
the asymptotic behavior of the perturbative series should set in and it cannot be 
ruled out that we are approaching this regime. Still, we would like to point out 
the smallness of this splitting compared to the total magnitude of the correction 
achieved by the reorganization of the perturbative series, which can be estimated 
by comparing the Coulomb line versus the NNNLO curve. In this respect, even 
if we consider the splitting between the NNLO and NNNLO curves as an error, 
its magnitude is rather small compared with the total gap. From this analysis, 
we conclude that we should use the RGI potential instead of the FO one and we 
will take this attitude in the rest of the paper. 

Another possibility that we explore is the use of the RS or RS' potential |iO] . 
To avoid numerical instabilities, due to the behavior of the potential at long 
distances, we also modify the potential in the following way: 



{VsD + 2(5mRs) I = 
{VsD + 25mRs)|^=i/r 



"+i(l/r) 



n=0 

oo 



n=0 



if r > Hr 
if r < fij. 



(60) 



Irrespectively of the potential we use, the short distance behavior of the potential 
and, consequently, yl('"^(ro; yu) is the same. The full expression for AW(ro;/i) IS 
more complicated in these cases than in Section 16.21 and we refrain from giving 
the general explicit expression and only show (for illustration) how it would look 
like at the lowest order. If for instance we consider LL running at short distance, 
namely 



r 1 



2tt 



■ In(yur) 



we have 



- 2mrCFas{fi)v{lo) 



27r 



/3oas(/i) 



/ 



IE + 



27T 



k 



f 



1E + 



27r 



(61) 



(62) 



(63) 



with f[x] = e^Ei(— s) — In x. The coefficients Vn have an expansion in and Iq 



v{lo) 




2(7s + l)/o + /o} 



+ (^) {(8 + 87.. + 47|)/o-(4 + 47^)/o^ + ^/o^} + 



(64) 



19 



-0.1 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 
H [GeV] 



2.0: 

1.5 : 




1 2 3 4 5 

ju [GeV] 



1400 ; 
1200 ; 




20 40 60 80 100 

ju [GeV] 



Figure 3: B^yjo) as a function of /i at LO (yellow), NLO (green), NNLO (blue) 
and NNNLO'(red) with fi, = 1.5 GeV, /i^ = 1 GeV and fip = fi^s = 0.7 GeV 
for charm, fig = fi^ = fip = GeV and fius = 1 GeV for bottom, and fig = 
Hr = Hf = '^^ GeV, and Hus = 10 GeV for top. Solid lines are obtained in the 
PS scheme using Eq. (l59l) and dashed lines are obtained in the RS' scheme using 
Eq. fl60|) . For reference we also include By^ (short-dashed black line). 
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We now perform the numerical evaluation of -B|^{o) at different orders in the 
static potential and compare the results obtained using the PS and RS' scheme. 
The results are shown in Figure [31 The difference between the schemes is small 
and converging for the case of bottom and top. In these two cases the differences 
between both schemes is pretty small for the NNNLO curves. This is again a 
good signal, since the dependence on the scheme is an indirect measure of the 
dependence on the long distance tail of the potential. For charm, the situation is 
less convincing. The gaps between schemes show marginal convergence at best as 
we increase the order. Yet, this gap is still much smaller than the gap between the 
Coulomb result and the NNNLO result. Comparing the Coulomb result, shown 
as the black short-dashed line, to our results, we can see that in all three cases 
a rather significant portion of the correction is already achieved with the LO 
RGI potential. In the case of top the NLO RGI potential is already quite close 
the most accurate NNNLO result. This behavior is also seen, to a lesser extent 
in the case of bottom. Note that the LO RGI potential exactly incorporates 
the r-dependent leading logarithms. This is equivalent to introducing an infinite 
number of corrections to the Coulomb potential and to iterate them an infinite 
numbers of times. This reorganization of perturbation theory seems to produce a 
major effect. Another observation is that the RS' scheme produces an accelerated 
convergence to the asymptotic regime. This is clearly seen in the top case, and to 
a lesser extent, in the bottom case. In those cases the LO RGI potential produces 
the bulk of the correction and the magnitude of the higher order corrections is 
smaller in the RS' than in the PS scheme. The price paid is that the splitting 
between different orders in the RS' scheme is less convergent. 

The dependence of the results on Hr is very small. Changing Hr from 2 GeV 
to 4 GeV for example results in differences that are an order of magnitude smaller 
than the changes we find by going from say LO to NLO. 

The dependence on fi will have to be cancelled by the scale dependence of 
the matching coefficient Cs{fi). Note that our evaluation of B^yju) also includes 
subleading logarithms, which are not matched by the precision of the RG (hard) 
computation. The fact that the scale dependence roughly corresponds to the 
Coulomb potential (with RG running) can be taken as an indication that sub- 
leading logarithms are not very important (see Figure [3] for illustration). 

Finally, there is also a dependence on the scale /i<j. This dependence (as the 
dependence on the renormalon subtraction scheme) partly reflects the dependence 
of the result on the long distance tail of the potential. On the other hand one 
can not take fig very small otherwise as(/is) becomes very large. We now perform 
the numerical evaluation of B^yio) at different orders in the loop expansion in the 
PS scheme and using different values of fig- The results for varying values of 
fig are shown as bands in Figure HI We also show the Coulomb result as the 
band enclosed by black dashed lines. This plot also illustrates that the bulk of 
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the correction is already achieved with the LO/NLO RGI potential in the top 
and bottom case, where we have convergence (in the charm case convergence 
is marginal at best). The fig dependence tends to diminish as one increases the 
order of the RGI potential in the top, and to a lesser extent in the bottom case. In 
the charm case the fig dependence remains almost constant. Overall, we find that 
the fig dependence is slightly larger than the scheme dependence, but still smaller 
than the typical gap due to working at different orders in the RGI potential. 

7 Phenomenology of the decay ratio 

Using the results obtained for B^lo) {E^^ ; /i) we can get improved determinations 
of the decay ratio, by combining Eqs. dHD, f|T^ and f H^ with the determination 
of Cg from Ref . [12] . We use the results obtained in Section 16.31 with the RGI 
potential, since they both achieve the resummation of logarithms and renormalon 
cancellation. The main source of uncertainties in the evaluation of -B^(o) {E^^ ; /i) 
is refiected by the computations at different orders in in the static potential 
and, to a lesser extent, by the dependence on iig. In comparison, the dependence 
on the quark mass, /ir, /i/ and ^us is small. Therefore, we will fix those parameters 
to the values used in Section 16. 3[ In Section 16.31 we also saw that the scheme 
dependence for renormalon subtraction was small, compared with the uncertainty 
due to the computation at different orders. Therefore, we will only take one 
scheme (PS) for reference in the plots. 

In order to explore different power counting expansions for our results, we 
will consider and compare different approximations. In particular we will show 
the effect of resumming logarithms in the matching coefficients D^-j g and Cg. 
We will see that the RGI in the matching coefficients plays an important role 
to make the result more factorization scale independent. The results obtained 
within a strict perturbative expansion (see Ref. pjij) are labelled as LO, NLO 
and NNLO respectively and, after resummation of logarithms, as LL, NLL and 
NNLL. Taking into account the static potential exactly, using numerical methods 
as described in the previous sections, we obtain improved predictions for the 
relativistic corrections that we label by including "I" to the previous labelling: 
NLLI (including Cg with NLL precision and the improved relativistic correction 
5pn) and NNLLI {cg with NNLL precision and the improved relativistic correction 
5pn)- For comparison we will also consider the result without resummation of the 
logarithms in the matching coefficient, NNLOI {cg with NNLO precision and 
the improved relativistic correction (5p„). For both, NNLLI and NNLOI we will 
consider the results taking the RGI static potential at LO, NLO, NNLO and 
NNNLO. 

From the point of view of a double counting in and v the NLL result (with 
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Figure 4: B^}o) using Eq. fl59|) as a function of /i at LO (yellow), NLO (green), 
NNLO (blue)' and NNNLO (red) with /i,. = 1 GeV and /i^ = fXus = 0.7 GeV for 
charm, /i^ = /xf = 2 GeV and fius = 1 GeV for bottom, and Hr = hf = '^^ GeV, 
and ^us = 10 GeV for top. The bands are obtained by variation of fig in the range 
1-1.5 GeV, 2-4 GeV and 20-60 GeV for charm, bottom and top respectively. For 
reference we also include By^ (grey band). 
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NLL precision for Cg) can be considered as 0{as,v'^) whereas NLLI is Olasjv"^) 
and NNLLI is 0(ag, f ^). As a general trend, moving from NLL to NLLI improves 
the scale dependence. This is due to the fact that, by using the RGI, NNLO 
0{a'^) logarithms count as NLL and can be matched with a part of the scale 
dependence of the relativistic 0{v'^) correction. Note as well that the inclusion 
of Cg with NNLL precision accounts for 0{a^) leading logarithms and beyond. 
Those should be cancelled by the inclusion of the subleading scale dependence 
of the relativistic correction. Most of it is actually built in by the numerical 
evaluation of the relativistic correction with the RG potential. In principle, this 
should be reflected in an improvement in the scale dependence in going from NLLI 
to NNLLI. On the other hand, this double counting in ctg and v scheme produces 
an unmatched scheme dependence, which can only be matched by working at the 
same order in ctg and v. 

We have also studied the dependence on the specific r-renormalization scheme 
of This depend ence should vanish when combined with cj^^. In 

particular we have studied the effect of eliminating 'Je in the logarithms in 
Eq. f l35|) and consequently using Eq. dH]) for c^^. Note that this is actually 
equivalent to using B^^lo^{E^^; jj,e~'^'^). We have checked that (at least in the 

cases where the series converges) this dependence fades away when considering 
the potential with increasing accuracy. The reason is that the je terms that 
appear at higher orders get more accurately described as we increase the order 
of our computation. This increases our confidence in the perturbative approach. 
The introduction of 'Je in the scale n makes the different terms in the expansion 
approach the asymptotic result faster, but the effect is not very significant. 

In the following subsections we will consider in turn the cases of top, bottom 
and charm. 

7.1 Top 

We start with the top since it is the cleanest possible case, where we expect best 
convergence. The scales are fixed as /ij? = /i^ = /i^ = 20 GeV and = 10 GeV 
and we work in the PS scheme. 

In Figure |5] we show the decay ratio at NNLOI (dashed lines) and NNLLI 
(solid lines) at different orders in ag in the static potential (LO: yellow; NLO: 
green; NNLO: blue; NNNLO: red). For reference we also include the LL, NLL, 
and NNLL results (short-dashed lines) obtained within a strict perturbative ex- 
pansion. Comparing the NNLOI with the NNLLI curves, it can be seen that the 
inclusion of the RG matching coefficients has a significant impact in reducing 
the scale dependence. Also, there is a sizable gap when moving from NNLL to 
NNLLI even if we take the LO RGI static potential (which includes the r run- 
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Figure 5: Decay ratio in the PS scheme at NNLOI (dashed) and NNLLI (sohd) at 
different orders in in the static potential {0{as): yellow; green; 0{a^): 

blue; 0(0;^): red). For reference we also include the LL, NLL, and NNLL results 
(short-dashed). 



ning producing the shift we observe in the plot). The inclusion of subleading 
corrections to the potential produces a convergent effect. Actually, the NLO RGl 
static potential result is already quite close to the asymptotic result. This may 
allow to define a counting in v, by taking the asymptotic limit of the series. The 
potential problem is that this counting in v is scheme dependent. 

To study this scheme dependence, in Figure M we show the decay ratio at 
NLLl in the MS and hard-matching scheme (see Ref. |12j and Eq. fHTj) ) and at 
NNLLI, all of them at 0{af) in the static potential. These results are compared 
to the LL, NNL and NNLL results. Moving from NLL to NLLl improves the scale 
dependence no matter what scheme is used. As we have already discussed, this 
is due to the fact that, by using the KG, NNLO O^a"^) logarithms count as NLL 
and can be matched with a part of the scale dependence of the relativistic 0{v'^) 
correction. On the other hand there is a sizable gap between the NLLl result 
obtained in the MS and hard-matching scheme. The latter is much closer to the 
full NNLLI result. The reason is that the two-loop hard correction is much smaller 
in the hard-matching scheme compared with the MS scheme. This could indicate 
that the hard-matching scheme leads to a more convergent series but it cannot 
be ruled out that this smallness is accidental for 0{a^). Therefore, we believe 
this gap gives a conservative estimate of the remaining uncertainties. Note that 
it is much larger than the other sources of uncertainties considered in this paper. 
For instance, we have also investigated the fig dependence and observed that it 
gets smaller when we consider higher orders in the static potential, pointing to 
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Figure 6: Decay ratio in the PS scheme at NLLI in the MS (grey dashed) and 
hard- matching scheme (grey sohd) and at NNLLI (red sohd). For reference we 
also include the LL, NLL, and NNLL results. 

the fact that the long-distance tail of the potential does not have a significant 
impact on the determination of the decay ratio. A similar comment applies to the 
renormalon scheme dependence. Therefore, in summary we find nice convergence 
in the top quark case. 

7.2 Bottom 

Turning to the bottom case, in Figure [7] we show the decay ratio in the PS 
scheme at NNLOI and NNLLI at different orders in as in the static potential. 
For reference we also include the LL, NLL, and NNLL results. We use fir = fJ'F = 
= 2 GeV and fius = 1 GeV. Again we can see that the inclusion of the RG 
matching coefficients has a significant impact in reducing the scale dependence. 
As in the top case, there is a sizable gap when moving from NNLL to NNLLI. 
The bulk of it is already obtained by taken the NLO(LO) RGI static potential 
in the PS(RS') scheme. The inclusion of subleading corrections to the potential 
produces a smaller effect, yet sizable. Compared to the top case the magnitude of 
the corrections is larger and the convergence using the static potential at different 
orders is worse, in particular in going from the 0(0^) to the Ola^) approximation 
of the static potential. Nevertheless, one can still see a band (though much wider 
than for top) where to roughly define a counting in v. We should also stress that 
using the 0{a^) RGI potential has some ambiguities, since ultrasoft effects enter 
at this order. Therefore, it can not be considered complete. 
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Figure 7: Decay ratio in the PS scheme at NNLOI (dashed) and NNLLI (sohd) at 
different orders in as in the static potential {0{as): yellow; 0{al): green; C(af): 
blue; (9(ag): red). For reference we also include the LL, NLL, and NNLL results 
(short-dashed). 

We study the scheme dependence in Figure [HI showing the decay ratio at 
NLLl in the MS and hard-matching scheme and at NNLLI. In all cases the static 
potential is taken at 0{a^). These results are compared with the LL, NNL and 
NNLL results. The general pattern of the results is similar to the top case. Mov- 
ing from NLL to NLLl improves the scale dependence irrespective of the scheme 
used. However, there is a sizable gap between the NLLl result obtained in the MS 
and hard-matching scheme, the latter is much closer to the full NNLLI result. As 
for top, we take this gap for a conservative estimate of the remaining uncertainty. 
Again, this gap is larger than other sources of uncertainties considered in this 
paper, like the splitting associated to different orders in the static potential, the 
lis or renormalon scheme dependence. Either way the errors are obviously larger 
here than in the top case. In particular we have found a larger sensitivity to [ig 
and the specific implementation of the initial conditions. 

We use this analysis to obtain an updated prediction for V{rji,(lS) — )■ 77). 
For the central value we use the NNLLI result with /i = 2 GeV and the set 
of parameters described before, obtaining 0.544 keV. The theoretical error has 
been estimated considering the difference between the NLLl (in the MS) and 
NNLLI result for /i = 2 GeV. We obtain 0.146 keV for this error. As we have 
already mentioned, we have checked that the uncertainties due the variation of 
these parameters, the scheme, or the consideration of different order in Us in 
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Figure 8: Decay ratio in the PS scheme at NLLI in the MS (grey dashed) and 
hard- matching scheme (grey sohd) and at NNLLI (red sohd). For reference we 
also include the LL, NLL, and NNLL results. 

the potential, is much smaller than the error quoted. Another source of error is 
experimental, coming from r(T(lS') e~^e~) = 1.340 ± 0.018 keV [H]. This 
produces a very small error: ±0.007 keV. Finally, we have also computed the 
error associated to the indetermination of as{Mz) = 0.118 ± 0.003. This error is 
even smaller: lo;oo4 keV. We combine the last two errors in quadrature and add 
linearly to the theoretical error (which completely dominates the error). After 
rounding we obtain T{rib{lS) — )■ 77) = 0.54 ± 0.15 keV. 

7.3 Charm 

Finally we consider the charmonium ground state. The applicability of our weak 
coupling approach to this system is doubtful. Nevertheless, we will find it reward- 
ing that the reorganization of the perturbative expansion significantly improves 
the agreement with the experimental data. Again we will use the PS scheme and 
set /is = 1.5 GeV, /ir = 1 GeV and fip = 1 = fius = 0.7 GeV. 

In Figure |9] we show the decay ratio at NNLOI and NNLLI at different orders 
in as in the static potential. For reference we also include the LL, NLL, and 
NNLL results. The experimental result, using r{J/ip — )• e'^e") = 5.55 ±0.14 keV 
and T{ric — ?■ 77) = 7.2 ± 0.7 ± 2.0 keV ^T] is shown as the light blue band, with 
the central value indicated by the horizontal solid light blue line. Once more 
we can see that the inclusion of the RG matching coefficients improves the scale 
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Figure 9: Decay ratio in the PS scheme at NNLOI (dashed) and NNLLI (sohd) 
at different orders in ag in the static potential (C(as): yellow; O^a"^): green; 
0{a^): blue; 0{a^): red). For reference we also include the LL, NLL, and NNLL 
results (short-dashed). The light blue band represents the experimental error of 
the ratio where the central value is given by the horizontal solid line. 

dependence and there is a sizable gap when moving from NNLL to NNLLI. The 
inclusion of subleading corrections to the potential produces a slightly smaller 
though still quite large effect. Compared to the bottom case the magnitude of 
the corrections is larger and the convergence is worse. We find the same problem 
in the associated evaluations of the energy and the wave function at the origin. 
Despite these shortcomings, the effect goes in the direction of bringing agreement 
with experiment. 

We study the scheme dependence by showing the decay ratio at NLLl in the 
MS and hard-matching scheme and at NNLLI in Figure [lOl The static potential 
is taken at 0{a^). The discussion is pretty similar to the top and bottom case. 
Moving from NLL to NLLl improves the scale dependence no matter what scheme 
is used. On the other hand there is a sizable gap between the NLLl result obtained 
in the MS and hard-matching scheme, the latter being much closer to the full 
NNLLI result. The reason is the same as for top and bottom. Taking this 
gap for an estimate of the typical size of the uncertainties produces an error of 
around 50% in the ratio. This encodes most of the experimental band and it is 
significantly larger than the typical spht produced by working at different orders 
in as in the static potential. 
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Figure 10: Decay ratio in the PS scheme at NLLI in the MS (grey dashed) and 
hard- matching scheme (grey sohd) and at NNLLI (red sohd). For reference we 
also include the LL, NLL, and NNLL results and the experimental ratio. 

8 Conclusions 

We have considered a different power counting in potential NRQCD by incor- 
porating the static potential exactly in the leading order Hamiltonian. In this 
scheme we compute the leading relativistic corrections to the inclusive electro- 
magnetic decay ratios. The effect of this new power counting is dramatic for 
charm, large for bottom, and sizable even for top. In the case of bottom, we 
produce an updated value for the rjb decay to two photons 

T{ribilS) 77) = 0.54 ± 0.15 keV. (65) 

In the case of charmonium, this scheme brings consistency between the weak 
coupling computation and the experimental value of the decay ratio, but the 
theoretical error is large. 

It is worth emphasizing that in the case where our expansion is more reliable, 
i.e. the top and bottom case, the bulk of the correction comes from using the 
first two orders of the RGI potential. The effect of higher-order corrections in the 
RGI potential is relatively small. The details of the importance of higher-order 
corrections depends on the scheme. In the RS' already the LO RGI potential 
gives the bulk of of the correction whereas in the PS two terms in the expansion 
are needed. Irrespectively, they both converge as one goes to higher orders. 

This approach could open the possibihty to reorganize the perturbative series 
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in a controlled way. We stress again that this is also relevant for top. Therefore, 
it is not a strong coupling effect but rather reflects the need of a more optimal 
resummation of perturbation theory. This might call for a reanalysis of previous 
results in this new scheme. It is an open question whether there is a similar effect 
in the case of the hyperfine splitting. We leave this discussion for a forthcoming 
paper. 

It would be misleading to only assign a theoretical error from the scale depen- 
dence. This is particularly obvious in the charm case, where the scale dependence 
by no means reflects a reasonable estimate of the size of higher-order corrections, 
which are difficult to estimate and can only be inferred from the apparent con- 
vergence of the expansion. 

Our formalism is flexible enough so that, with little effort, we could replace the 
perturbative static potential by any potential, in particular, by one fitted to non- 
perturbative lattice data. This could be of particular relevance for charmonium 
but it could also be of help for bottomonium, provided the static potential is know 
with enough accuracy in the unquenched approximation. This would eliminate 
the error associated to higher order terms in the static potential, but not the error 
due to higher order terms in the hard matching coefficient and the associated RG 
improvement. 
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